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We study the quantum phase transition mechanisms that arise in the Interacting Boson Model. 
We show that the second-order nature of the phase transition from U(5) to 0(6) may be attributed 
to quantum integrability, whereas all the first-order phase transitions of the model are due to level 
repulsion with one singular point of level crossing. We propose a model Hamiltonian with a true 
first-order phase transition for finite systems due to level crossings. 
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Quantum phase transitions (QPTs) have attracted 
great attention from the theoretical and experimental 
communities in recent years. Experiments on high- Tc 
superconductors, on quantum-Hall systems, on Bose or 
Fermi dilute gases in different trap geometries and on 
nuclei close to a critical point challenge theoreticians to 
develop reliable approaches with which to describe the 
critical properties of these systems. 

A QPT describes a structural change in the properties 
of the ground state of the system associated with the vari- 
ation of a control parameter (the temperature being the 
usual one in classical phase transitions) . This parameter 
might be the level of doping in high- Tc superconductors, 
the magnetic field in quantum-Hall systems, the scatter- 
ing length in dilute gases, the nucleon number in nuclear 
phase transitions, etc... 

The phase diagram of the Interacting Boson Model 
(IBM) [1] of nuclei has been studied at the mean-field 
level [2,3], using both Catastrophe theory [4] and the 
Landau theory of phase transitions [5] . Both treatments 
show that the model displays first-order phase transi- 
tions from spherical to deformed shapes and from oblate 
to prolate deformed shapes, with an exceptional point of 
an isolated second-order phase transition. In this work, 
we study in detail the phase transitions associated with 
the IBM. 

The IBM is a phenomenological model of nuclear struc- 
ture with a deep connection to the underlying micro- 
scopic Shell Model. In this model, which has been very 
successful in describing nuclear spectroscopic properties 
of collective nuclei, pairs of correlated nucleons of angu- 
lar momentum L — and L = 2 are represented by s and 
d bosons, respectively. In its simplest version the model 
does not distinguish between protons and neutrons and 
has an underlying group structure based on the dynami- 
cal group U(6). 

The most general IBM Hamiltonian can be written in 
terms of 6 free parameters, as a linear combination of the 
linear and quadratic Casimir operators of U{5) and the 
quadratic Casimir operators of 0(6), SU{3), 0{5) and 
0(3). A convenient and frequently-used form of the IBM 
Hamiltonian that keeps all the main ingredients related 



to the structure of the ground state is the Consistent-Q 
Hamiltonian [6] 
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N is the total number of bosons, which is equal to the 
number of nucleon pairs in the valence space. 

One of the most important features of the IBM is the 
existence of four distinct dynamical symmetries (DS), 
each representing a well-defined phase of nuclear collec- 
tive motion. A quantum system has a DS if the Hamil- 
tonian can be expressed as a linear combination of the 
Casimir operators of a subgroup chain of the dynamical 
group. The four dynamical symmetries are: the U{b) 
symmetry for spherical vibrational nuclei [x = 1), the 
S'J7(3) symmetry for prolate deformed nuclei [x = 0, 
X — — \/7 /2), the 0(6) symmetry for 7 unstable deformed 
nuclei (x = 0, x = 0) and the SU (3) symmetry for oblate 
deformed nuclei (x = 0, x = v^/2). 

The symmetries of the IBM and the transitions be- 
tween them are illustrated in the extended Casten trian- 
gle [7] of Figure 1, along with the corresponding values of 
the parameters x and x- Distinct phases associated with 
three "shapes" , denoted S for spherical, P for prolate de- 
formed and O for oblate deformed, can be seen in the 
figure. The three different shape phases are separated 
by two lines of first-order phase transitions, a; = 0.8 and 
X = 0, respectively, with the exception of a triple point 
[5] at a; = 0.8 and x = O7 which is represented by a soHd 
square and for which the phase transition is second or- 
der. Recently, two effective models based on the Bohr 
Hamiltonian, called A(5) and i?(5), have been proposed 
to describe the physics at the critical points along the 
C/(5) - SU{3) and U{b) - 0(6) lines, respectively [8]. 

The aim of this letter is to go beyond the mean-field 
description of the phase diagram in Fig. 1. We will show 
that the first-order QPTs arise due to level repulsion be- 
tween two states with different properties, as was earlier 
recognized in [9]. However, there are two exceptional sin- 
gular critical points, represented by the bold square and 
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FIG. 1. Phase diagram of the Interacting Boson Model in 
the phase space of the control parameters x and x- S, O, and 
P stand for the spherical, oblate, and prolate phases respec- 
tively. 

the bold diamond in Fig. 1, for which the mechanism 
triggering the phase transition is completely different. 

We begin with a discussion of the transition from 0(6) 
to U{5). It is described within the Hamiltonian (1) by 
setting X = varying x from to 1. Earlier investi- 
gations of the properties of the level spacings along this 
leg of the Casten triangle showed a Poisson distribution. 
This indicates that the model is intcgrablc along this leg 
[10], even though it is not characterized by a global Dy- 
namical Symmetry. 

Quantum integrability requires the existence of a com- 
plete set of mutually commuting hermitian operators. 
These operators are the integrals of motion or quantum 
invariants and their eigenvalues are the conserved quan- 
tities that completely label a unique basis of common 
eigenstates. A DS is a particular class of integrable model 
in which each one of the Casimir operators in the sub- 
group decomposition chain is a quantum invariant. But 
there are also other integrable and exactly solvable mod- 
els that have played important roles in the understanding 
of quantum many-body systems. Examples include the 
Heiscnberg and Hubbard models, both of which are solv- 
able by the Bethe ansatz. The Hamiltonian (1) for x = 
is a particular case of a general class of pairing models 
that have been shown to be quantum integrable [11]. The 
quantum invariants for boson pairing Hamiltonians and 
the corresponding exact solutions were discussed in Ref. 
[12]. Here we will specify the form of the quantum in- 
variants in the 0(6) — U{5) transition. Since we are only 
interested in properties of the ground state, we focus our 
treatment on the case in which all bosons are paired to 
zero angular momentum. 

We first note that this transition is embedded in the 
decomposition U{6) D ■■■ D 0(5) D 0(3) D 0(2), where 
a missing parameter dependent "symmetry" or quantum 
invariant commuting with the Casimir operators of all the 



FIG. 2. Transition from 0(6) to U{5). The ground state 
energy (left panel), expectation value of the nd (central 
panel), and derivative of the ground state energy (right panel) 
are presented as a function of the control parameter x for a 
system with 1000 bosons. The exact results are represented 
by circles and the mean field results by a full line. 

groups in the chain is responsible for the integrability of 
the model. 

As noted in references [10,13], the Hamiltonian itself, 
interpolating between the Casimir operators of f/(5) and 
0(6), constitutes the additional independent quantum 

invariant ensuring the integrability of the system. 

The eigenvalues of the Hamiltonian (1) for x = are 
[12] 

a— 1 

where the pair energies ea are the solutions of the coupled 
set of nonlinear Richardson equations 



2(1 -x) e„ 2-ea e.a-e.0 

The Hamiltonian eigenvalues (2) are analytic functions 
of the control parameter x. Therefore, in the absence 
of level crossings at the ground state level, a first-order 
phase transition is precluded. Level crossings can occur 
only for integrable models with more than one parameter 
dependent integral of motion, an example being the sdg 
boson model [14]. 

Because of the exact solvability of the model, we are 
able to obtain exact results numerically for a large num- 
ber of bosons. In Fig. 2, we show results for 1000 bosons 
and compare them with the mean-field solution. In the 
left panel, we plot the ground-state energy versus the 
control parameter x; in the central panel, the expecta- 
tion value of the number of d bosons < > is presented 
as a function of x\ in the right panel, we plot the deriva- 
tive of the energy versus x. Due to the exact solvability 
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FIG. 3. Ten lowest 0^ states in the transition from U{5) 
to SU{3) for N = 50 bosons as a function of x. On the right 
panels we show in an enlarged scale two examples of level 
repulsion. 



of the model, we can make use of the Hellmann-Feynman 
theorem to express the first derivative of the energy as 
the expectation value of the derivative of the Hamiltonian 
(1) with respect to x, i.e. || = + j^Q° ■ Q°. Both 
operators, and ^ , are zero in one phase (disordered) 
and different from zero in the other phase (ordered), ful- 
filling the requirements for being order parameters. Both 
expectation values are continuous at the transition point 
satisfying both the Landau and Ehrenfest definitions of 
a second-order phase transition. 

We will now discuss the phenomenon of level repulsion 
for non-integrable systems. For an arbitrary finite value 
of X, we select a value of x = a;o in the vicinity of which 
we can use non-degenerate perturbation theory. Defin- 
ing x' = X — .xq, we can rewrite the Hamiltonian as H = 
xorid + ^Q"" -Q^ + x' {nd + j^Q^ ■ Q^) =Hq+ x'V . 
The first and second derivatives of the Hamiltonian eigen- 
values in perturbation theory are [15] : 



element V12 remains finite as A ^ 0, there is an infinite 
repulsion and this prevents the particles from colliding. 
This is a pictorial view of the phenomenon of level repul- 
sion for non-integrable Hamiltonians. On the other hand, 
if there is an additional conserved quantity and both lev- 
els have different symmetries, then V12 is strictly zero 
and the two levels can cross. In what follows, we discuss 
the application of these ideas first to the ?7(5) to SU (3) 
transition and then to the SU{?>) to SU{?>) transition. 

In Fig. 3, we show in the left panel the ten lowest 
0+ states for a system of 50 bosons as a function of x 
along the C/(5) to 5J7(3) transition. There are several 
cases in which two 0+ levels get very close in energy but, 
as just discussed, the closer they get the larger the level 
repulsion between them is. In the right panels, we mag- 
nify two of the cases of closest approach to show that 
indeed the levels repel and do not cross. In particular, 
there is a smooth crossover in the ground-state energy at 
the point of closest approach with the first excited state. 
The dynamics of these two levels is well described by Eq. 
(6) with finite matrix elements Vij. The crossover will 
change gradually as the number of bosons increases to- 
wards a non-analytic point in the thermodynamic limit, 
defining a first-order phase transition. From the above, 
we conclude that level repulsion is indeed responsible for 
the first-order phase transition in this case [9] . 

Now we turn to the third leg of the triangle (a; = 0), 
the transition from deformed oblate {SU{3) shapes, x = 
V7/2) to deformed prolate {SU{3) shapes, x = -V7/2) 
shapes, passing through the 7-unstable 0(6) case with 
X = 0. In this case, the first-order phase transition takes 
place precisely at the 0(6) dynamical symmetry [7]. At 
this critical point, the 0+ states can be classified with the 
T quantum number of an 0(5) symmetry. As discussed 
above, the realization of this new symmetry allows the 
crossing of levels with different quantum numbers. The 
Hamiltonian (1) at this point is 
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(C2[0(6)]-C2[0(5)]) , (7) 
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Assuming that there is a pair of energy levels that 
get very close in energy, the energy gap between them 

is A {x) = E2 (x) — El (x) . Neglecting the contribution 
of all other states, perturbation theory for the gap deriva- 
tive gives 



OX 



Vn, 



dx^ 



(6) 



We can interpret Ei and E2 as the positions of two 
fictitious classical particles in a one-dimensional space. 
Then (6) are the equations for their relative velocity and 
the force between them. Clearly, if the interaction matrix 



with eigenvalues E^yr = [-o" (cr + 4) + r (t + 3)] /2N. 
The lowest 0(6) band corresponds to <j = TV and r = 
0, 3, • • •. In the large N limit, the 0(5) states within the 
ground state 0(6) band will be degenerate, since the en- 
ergy spacings go as N~^. According to the above anal- 
ysis, in the thermodynamic limit an infinite set of 0"*" 
states with different seniority quantum numbers r cross. 
The crossing is allowed because of the 0(5) symmetry at 
the critical point, implying that the interaction matrix 
elements in (6) are zero. 

As just noted, the exact degeneracy in the levels of the 
ground-state band at .t = and X = only occurs in 
the large-N or thermodynamic limit. To better visualize 
this QPT in the context of finite systems, we therefore 
consider a hamiltonian similar to that in (1), but with an 
additional term that cancels the 0(5) term in it that pro- 
duces energy spacings within the ground-state band for 
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FIG. 4. Transition from SU{3,) to .S'r/(3) through the 0(6) 
limit. All 0^ cigcnstatcs of the Hamiltouian (8) for a system 
with N = 9 bosons are represented as a function of the control 
parameter x- The right panels show in a larger scale the level 
repulsion of two close levels at x (right-top panel) and 
the level crossing at x = (right-bottom panel). 



finite N (see eq. (7)). With this in mind, we consider the 
following Hamiltonian associated with the a; = transi- 
tion leg: 



H : 



1 

'n 



V7 



C2 [0(5)] 



(8) 



For X = the Hamiltonian (8) is in the SU{i) 

limit; for x — it is in the SU{i) limit; and for x = 
it is proportional to the quadratic Casimir operator 
of the 0(6) group. The second term in the modified 
Hamiltonian (8) thus cancels out the 0(5) breaking of 
the degeneracies in the ground-state band that occurs 
for finite N. 

In Fig. 4, we show all of the 0+ eigenstates of the 
Hamiltonian (8) for a system with N = 9 bosons as a 
function of x- At x = 0, the four levels that correspond 
to the ground-state band with a = 9 cross. This crossing 
gives rise to a first-order QPT for a "finite" system. This 
is the unique point along the a; = log at which level 
crossings are allowed due to the realization of the 0(5) 
symmetry. Away from this point the mechanism of level 
repulsion sets in. In the upper right panel of Fig. 4, we 
show in an enlarged scale an apparent level crossing at 
X 7^ 0. The mechanism of level repulsion predicted in eq. 
(6) is numerically verified. On the contrary at x = 0, the 
realization of the 0(5) symmetry allows the crossing of 
levels with different r quantum numbers, as shown in the 
lower right panel. 

In summary, there are two singular points of QPTs in 
the IBM phase diagram. The first is in the transition 
from U{5) to 0(6) and corresponds to a second-order 
QPT due to the integrable nature of the Hamiltonian 
along this line, which ensures that the ground-state en- 
ergy is an analytic function of the control parameter. The 
second singular point is located at the QPT from SU{3) 



to SU{3), precisely at the 0(6) critical symmetry. In the 
thermodynamic limit, this QPT is associated with the 
crossing of an infinite number of 0+ states. By adding 
to the Hamiltonian (7) a term proportional to the 0(5) 
Casimir operator, we were able to realize this QPT for a 
finite system. This is an unusual, but especially interest- 
ing, case to study further, since most phase transitions 
take place in the thermodynamic limit. Therefore, the 
model described by the Hamiltonian (8) may be useful 
to explore the universal properties of first-order phase 
transitions. 

These two isolated critical points represent exceptional 
points in the phase diagram of the IBM. Along the lines 
separating the spherical, prolate and oblate phases, the 
mechanism underlying the first-order phase transitions 
is level repulsion between the ground state and the first 
excited state. For finite N, this mechanism is reflected 
by a soft crossover, but in the large N limit the model 
develops a singularity in the derivative of the gap consis- 
tent with a first-order phase transition. This conclusion 
is also valid for more general IBM hamiltonians, since 
level repulsion is a universal phenomenon associated with 
quantum non-integrability. 
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